Boundary-locality and perturbative structure of entanglement spectra in gapped systems 
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The entanglement between two parts of a many-body system can be characterized in detail by the entangle- 
ment spectrum. Focusing on gapped phases of several one-dimensional systems, we show how this spectrum is 
dominated by contributions from the boundary between the parts. This contradicts the view of an "entanglement 
Hamiltonian" as a bulk entity. The boundary-local nature of the entanglement spectrum is clarified through its 
hierarchical level structure, through the combination of two single-boundary spectra to form a two-boundary 
spectrum, and finally through consideration of dominant eigenfunctions of the entanglement Hamiltonian. We 
show consequences of boundary-locality for perturbative calculations of the entanglement spectrum. 



Introduction — The study of entanglement-related quan- 
tities in condensed matter has led to a large body of interdis- 
ciplinary work [1]. Recently, the concept of the entanglement 
spectrum has established itself at the forefront of the field, as 
this spectrum provides much finer information than a single 
number like the von Neumann entropy. Considering a biparti- 
tion of the system into parts A and B, the entanglement spec- 
trum (ES), {£,i}, is defined in terms of the Schmidt decompo- 
sition 
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Here | ■0 ) is the ground state, and the states | ipf ) ( | ipf ) ) form 
an orthonormal basis for the subsystem A (B). The ES {£,i = 
— log Xi} can also be thought of in terms of the eigenvalues 
{Xi} of the reduced density matrix pA obtained after tracing 
out the B part of the system density matrix \ip){il;\. 

The ES has been studied earlier for insights into the density 
matrix renormalization group (DMRG) algorithm [2-4], and 
more recently because of its relation to low-energy boundary- 
related modes in topological phases [5-7]. One way of think- 
ing of the ES is in terms of an entanglement Hamiltonian (EH) 
Ha = — log PA acting only on the A degrees of freedom 
[3, 9-12]. The ES is then the spectrum of this object. This 
gives a suggestive correspondence to the energy spectrum of 
true bulk Hamiltonians. This idea has also inspired calcula- 
tions of bulk terms that appear in such a 'Hamiltonian' [3, 9]. 
In topologically ordered systems such as fractional quantum 
Hall states, the EH indeed seems to share the low lying edge 
modes with the physical Hamiltonian of a system with the 
same boundary. This correspondence gave rise to the loose 
idea that entanglement Hamiltonians are quite similar to phys- 
ical "bulk" Hamiltonians, apart perhaps from non-essential 
renormalizations, like longer range hoppings or spatial inho- 
mogeneities [5, 6]. However, the well-known area law for the 
entanglement entropy Sa — ^'Yl ^i log -^i' stating that Sa 
scales with the size of the boundary between A and B parts, 
clearly suggests that the spectrum of pA (and hence also the 
ES and EH) must in some sense be dominated by the boundary 
degrees of freedom. 

This contradiction is highlighted even more by considering 
states that are described by simple matrix product states, such 



as the Majumdar-Ghosh point in the frustrated antiferromag- 
netic chain (J1-J2 model) or the AKLT (Affleck-Kennedy- 
Lieb-Tasaki) state [13]. The entanglement spectrum for these 
states consists only of small number of finite values, inde- 
pendent of the (sufficiently large) block size A. These simple 
gapped states thus clearly have entanglement spectra deter- 
mined completely by the boundary and not by the bulk. 

In this work, we sharpen the boundary picture of the entan- 
glement spectrum by considering more general gapped one- 
dimensional (ID) systems. The boundary-local nature of the 
ES is demonstrated and made quantitative in several ways. We 
show that the ES of gapped states generally have a hierarchi- 
cal structure, with excitations farther from the boundary be- 
ing successively higher in the ES. This also leads us to the 
idea that the ES can be obtained from a boundary-connected 
peiturbative calculation, and we show elements of such a cal- 
culation. We also show that the ES of two-boundary blocks 
can be constructed by combining the ES of single-boundary 
blocks, demonstrating that the bulk degrees of freedom play 
a secondary role. We will concentrate on the widely famil- 
iar XXZ chain, and as a bonus, we uncover a beautiful set 
of degeneracy structures in the ES of this model. We believe 
these findings are very general; we emphasize this by briefly 
presenting boundary locality in two other systems, namely the 
Heisenberg ladder and the Mott insulating phase of the Bose- 
Hubbard model. 

In our examples we consider perturbations around a product 
state; it should be possible to adapt to perturbations around a 
point which is a matrix product state with small-rank matrices. 
Gapped ID phases often contain or are connected to such sim- 
ple points. Gapped phases in general are not very universal in 
their low-energy properties; it is therefore remarkable that the 
ES has features common in many gapped phases, which we 
identify in this work. 

The usual methods for calculating ES are numerical exact 
diagonalization or DMRG, or using two-point correlators for 
non-interacting systems [2]. We propose our boundary-linked 
peiturbation theory as a general alternate technique for calcu- 
lating ES in gapped phases containing a simple point around 
which one can perturb. The basic insight is that, while cal- 
culating the complete reduced density matrix peiturbatively 
requires perturbation terms acting on the whole system/block. 
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FIG. 1. (Color online) Single boundary entanglement spectrum (ES) of the XXZ chain, (a) Long or infinite blocks. The ES levels appear at 
well-defined levels, which correspond to successive perturbative orders n. The accompanying integers denote degeneracies. On the right we 
show the total degeneracies at each order. The levels marked (i) and (ii) have dominant contributions shown in (b), where a fiUed-rectangle 
notation is introduced for domain walls, (c) Dominant contributions to low-lying levels, (d) Dominant contributions to low-lying 5SX — 
levels, (e) Finite-size block (12 spins in A); arbitrary-precision exact diagonalization data for A = 10. The effects of finite size are seen 
through breaking of degeneracies at order ~ 8 and higher. The integer-indicated degeneracies are now approximate. 



if one is only interested in determining the ES and the lead- 
ing order eigenfunctions (entanglement eigenstates), then a 
boundary local perturbation theory is able to construct the ES 
order by order in a physically transparent way. 

Boundary-linked perturbation theory — To obtain ES lev- 
els correctly up to n-th order, one only needs to consider per- 
turbations up to this order which act within an sites from 
the boundary, a depending on the type of perturbation. This 
is a direct but powerful consequence of boundary-locality, 
namely, that the ES levels at increasing orders correspond to 
increasing distances from the boundary. The perturbative cal- 
culation necessary for calculating ES levels is thus "boundary- 
linked". 

For any two orthonormal bases \'ff')^ Ivf )' of subsystems 
A and B, the system wavefunction can be written as |-0) = 
Y, M,j |(^f ) (8) |(/3f ). The e"'^'/2 ^f^^ (1) are obtained from a 

id 

singular value decomposition (SVD) of the matrix M. Further 
aspects of the perturbative calculation of ES can be formulated 
in terms of the matrix M. For example, when calculating M 
perturbatively, we find that a contribution to M at some or- 
der can lead to a new ES level only if the contribution is not 
appearing on the same row or the same column as a previous 
contribution which led to a new ES level at lower order In 
physical terms, this means that we do not get new ES levels 
by applying a perturbation only to the A block but need to ap- 
ply perturbation terms to both blocks in a linked way. Some 
perturbative calculations for the Heisenberg chain, and some 
general perturbative features of the matrix M, appear in the 
Supplementary Materials. 

The XXZ gapped phase — We consider the anisotropic 
Heisenberg (XXZ) chain, H = Hz + T~Lxy, with T-Lz = 



AE^f^f+i. and H.y ^ ^ E(5,+ ^r+i + h.c), in the 

i i 

gapped phase A > 1. The Ising limit A ^ 1 is simple: 
the ground state is spanned by the two product (Neel) states 
|N1) = I titltl • • • ) and |N2) = | iUtlt • • • >■ We will 
therefore think perturbatively around this simple limit, T-L^y 
being the perturbation. In the language of 'domain walls' 
(bond between neighboring aligned spins), the perturbation 
T-Lxy can have two effects: it can create a pair of domain walls 
situated two bonds apart, or it can move a domain wall over 
two sites. We choose for simplicity to work with perturba- 
tions around a single Neel state, by selecting |N1 ) instead of 
the linear combination ^(|N1) + |N2)) to be the vacuum 
state. For cases where the linear combination is appropriate, 
the ES can readily be reconstructed from knowledge about the 
individual Neel states by overlaying and shifting by log 2 of 
the single state ES. 

The XXZ single-boundary ES — Fig. 1(a) shows the ES 
for the simplest setup, namely, a single boundary partitioning 
an open XXZ chain into two long blocks. The total S^' of the 
block A is a good quantum number for the reduced density 
matrix p^. Hence we organize the ES into sectors of 5S\, 
which is the difference of the S\ value of the corresponding 
entanglement eigenstate from the block 5*^ value in the Neel 
state (0 or ±i). The most notable feature is that the ES is 
organized in equally spaced hierarchical levels. In Figs. 1(a) 
and 1(e) the hierarchical orders are denoted as "perturbation 
order" because successive orders in a perturbative calculation 
yield successively higher ES levels. 

Higher levels correspond to excitations farther from the 
boundary. In Figs. 1(c) and 1(d) this is illustrated pictorially 
by showing the dominant configurations in some of the states 



occurring at several levels. 

The ES can be understood physically, and constructed ac- 
curately and completely, through the boundary-linked pertur- 
bation theory introduced previously. Here we present a qual- 
itative discussion highlighting the algebraic structure of the 
problem. For the XXZ model, this theory can be formulated 
as follows. The vacuum state has no domain walls. Excita- 
tions in the ES are obtained by introducing domain walls near 
the boundary and by moving them into the bulk. From ex- 
plicit perturbative calculations at the first few orders, we find 
that the only perturbation terms leading to new ES levels are 
those which do one of the following: (1) create a pair of do- 
main walls by applying Hxy across the boundary: 



(a) ...mt][wt. 
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(2) 



or (2) act twice, on the first bonds on either side of the bound- 
ary, again creating a domain wall pair: 



.am][m<^- 



.aiimn^- 



(3) 



or (3) act twice symmetrically with respect to the boundary so 
as to move a domain wall on each block, each by two bonds 
away from the boundary. Above, domain wall positions are 
shown underlined, and a represents a spin of unspecified ori- 
entation. The first process changes S^ by ±1, while the others 
preserve S%. The boundary-locality of the ES is encoded in 
the restriction that domain walls are created only through the 
first two processes above and not farther away from the bound- 
ary. The symmetrical pairwise application of Hxy ensures that 
the new contribution is to an element of the M matrix that is 
not on the same row or column as a previous element that al- 
ready led to a lower order ES level. 

Given a domain wall configuration, the corresponding A 
state belongs to the sector 6S\ = (-1)^°^ I]j(-1)-'(1 ^ 
(—1)^^), where N^ is the number of domain walls in A, 
j e [IjiV^i] labels those starting from the boundary, and 
£j > 1 is the position of the j-th domain wall measured from 
the cut. In our perturbative rules, the order at which a partic- 
ular configuration appears is given by n = ^ £j. One can 
thus construct the ES, organized by SS^ sector, from the per- 
turbative rules in terms of the domain walls. In this way one 
obtains exactly the ES of Fig. 1(a). 

For each SSj^, the lowest ES level (filled symbols in 
Fig. 1(a)) is generated by a state where the region near the 
boundary is packed with domain walls at each bond, creat- 
ing a ferromagnetic region as large as needed to have that 



6S% value. 



The lowest level for sector SS^ occurs at order 



m{SSl) = \5Sli2SS^A-l)\. 

The degeneracies of the ES follow intriguing patterns, 
which can be explained through our construction. The degen- 
eracy at order n in a given dS% tower is p{[n — m{SS%)]/2) 
where p{x) is the number of integer partitions of the integer 
X. The total degeneracy at a given order n, listed to the right 
of Fig. 1(a), is given by q{n), the number of partitions into 
unequal summands [14]. The degeneracy sequence q{n) was 
observed in corner transfer matrix calculations [15], but our 
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FIG. 2. (Color online) Combining two single-boundary ES to form 
a two-boundary ES, for the XXZ chain. Shading indicates the A 
partition. (a,b) The single-boundary ES obtained by considering two 
different bipartitions. (c) Combination of single-boundary cases to 
form a two-boundary ES (large empty diamonds), compared with 
numerical diagonalization data at A = 10 (small filled diamonds), 
which match perfectly. 



perturbative construction now gives a physical picture: the 
number of ES levels at order n is the number of ways one can 
place domain walls within the first ri positions while keeping 
the sum of position labels (^ £j) to be n; this number is q{n) 
by definition. 

For the infinite XXZ chain, this ES is actually exact for 
all A > 1 (not just A ^ 1), with the ES level spacing 
2arccoshA, as a result of integrability [15]. In a generic (non- 
integrable) gapped phase, the hierarchical level structure is 
most pronounced in the case of weak perburbations, and gets 
progressively broadened by higher order renormalizations of 
entanglement levels upon increasing the perturbation. 

Fig. 1(e) plots the exact ES for a finite open chain. The ES 
follows the infinite-system structure at low orders, and only 
deviates significantly (e.g. split degeneracies) at orders corre- 
sponding to the distance between the boundary and the phys- 
ical edge. This illustrates again that the lowest ES levels cor- 
respond to smallest distances from the partition boundary. 

The two-boundary ES — Boundary-locality implies that 
the single-boundary ES is sufficient to construct the ES 
for multiple-boundary blocks when the boundaries are suffi- 
ciently far apart. Formally, this physical intuition means that 
the reduced density matrix pA of a two-boundary block should 
'factorize' in the loose sense that p^ is iso-spectral to pl®Pr, 
where Pli^r) are (virtual) reduced density matrices appropri- 
ate for single-boundary blocks having the left (right) boundary 
of A. This idea has recently be shown to be at work for ES of 
blocks of a fractional quantum Hall state on the torus [6] . 

Fig. 2 shows the combination of two single-boundary ES to 
form a two-boundary ES in an XXZ chain. The case shown 
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FIG. 3. (Color online) Single-boundary ES (DMRG data) for the (a) 
Heisenberg spin ladder and (b) Bose-Hubbard chain. For the ladder, 
we show numerical data for Ji,.g = 0-lJ,u?ig, in the rung singlet 
phase. The partitioning is shown in the top cartoon (region A shaded 
as in Fig. 1). The integers next to the data points indicate approxi- 
mate degeneracies. In the associated cartoons, we show with shading 
the rungs (degrees of freedom) which contribute to various levels of 
the ES. For the Bose-Hubbard case (data shown for U — 40), the 
cartoons show dominant contributions to ES levels. Sites with oc- 
cupancy 7^ 1 are shaded; occupancies are indicated with integers on 
each site. Partition A is again as in Fig. 1 . In both panels the horizon- 
tal light dashed lines are guides to the eye indicating the perturbation 
order at which the levels appear. 



corresponds to the ground state being a dressed single Neel 
state |A^1), and an even size for the A block. The relevant 
single edge ES, Figs. 2(a) and 2(b), are then (5S'^)-inverted 
versions of each other Each level of Fig. 2(a) is combined 
with each level of Fig. 2(b), the two constituents contributing 
additively to the order and to SS^. The resulting ES has the 
structure and degeneracies shown in Fig. 2(c). In the com- 
bined ES, the lowest level at any dS% is found, from a combi- 
nation ofmL{x) — \x{2x + 1)1 and mji{x) = \x{2x — 1)|, to 
be m{SS^) — {6S^y. This explains the parabolic envelope 
shown as dashed line in Fig. 2(c). 

If the A block is odd instead of even, the two single- 
boundary ES are identical rather than reflections of each other; 
the combined ES would then not be symmetric under SS^ in- 
version. 

Other gapped ID systems — Boundary-locality can be il- 
lustrated with the ES of various other gapped systems. As a 
first additional example we demonstrate in Fig. 3(a) the case 
of the gapped spin-i Heisenberg ladder [16] with a cut which 
bisects both chains (in contrast to the chain - chain setup stud- 
ied in Refs. [10, 12, 18]). There is once again a hierarchical 
level structure corresponding to increasing distances from the 
boundary. The levels have SU{2) structure; there is a triplet at 
first order, no new levels at second order, two triplets at third 
order, a singlet and a quintuplet at fourth order, etc. Boundary- 
linked perturbation arguments can explain this structure, and 
also the parabolic envelope. 

Our last example system is the Bose-Hubbard chain in the 
gapped Mott insulating phase at unit filling (Fig. 3b). The 
level structures in the ES are richer than the XXZ or ladder 



cases, due to the larger local Hilbert space. Some features, 
like a near-parabolic envelope, appeared in Ref. [17]. Here we 
provide the physical picture behind this observation based on 
boundary-locality, as illustrated through the cartoons in Fig. 3. 
The nature of the local Hilbeit spaces now allow asymmetric 
excitations in the two blocks, leading to a modified parabola 
for the envelope: ^ fv x{x + l)\nU -ln[{x + 1)1]+0{1/U), 
with a; = \6Na\- 

Conclusions — While entanglement as a boundary- 
related quantity is a direct reflection of the area law and is thus 
well-appreciated, the competing picture of a 'bulk' entangle- 
ment Hamiltonian clearly calls for a thorough exploration of 
boundary-locality and its consequences. We have provided 
the elements of such an analysis here. In particular, we have 
shown how boundary-locality allows a perturbative formal- 
ism specific to calculating entanglement spectra. At any finite 
perturbative order, we reproduce a finite number of the lowest 
ES levels, which amounts to constructing a finite-rank matrix 
product representation of the state. 

To put our work in context of recent literature: high tem- 
perature series expansions for the mutual information exhibit 
a boundary linked property [19] similar to our rather different 
T = peiturbation theory targeting ES. Other recent works 
highlighted the structure of ES along a ID block boundary in 
2D systems [8, 18]. Our work is complementary in that we 
clarify the structure of the ES perpendicular to the cut, i.e. as 
one moves further into the bulk. 

Our work opens up several research avenues. For example, 
in gapless ID states, and in Fermi liquids in higher dimen- 
sions, the entanglement entropy does not scale purely with 
the boundary size but has a logarithmic correction. This im- 
plies some kind of weakening of the boundary-local picture; 
degrees of freedom deep in the bulk must play a greater role 
in such cases. Quantifying this effect remains an open task. 
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SUPPLEMENTARY MATERIALS 

In these Supplements, we 

(1) present general features of the perturbative structure of the 
'Schmidt' (M) matrix and the entanglement spectrum (ES), 
and a relatively general proof of boundary-locality at low per- 
turbation orders; 

(2) provide some details concerning the perturbative calcula- 
tion of ES in the Heisenberg spin ladder. 



OVERVIEW; THE 'SCHMIDT' (M) MATRIX 

The purpose of these Supplementary Materials is to pro- 
vide some details and examples of perturbative calculations 
of the entanglement spectrum (ES), both in specific and gen- 
eral settings. The general formulation also provides a proof of 
boundary-locality for the ES, for a wide class of gapped ID 
systems described by local Hamiltonians. 

It is convenient to describe our calculation in terms of the 
matrix ]V[, appearing in the Schmidt decomposition process, 
as introduced in the main text: 



|^) = ^M 



^JlVf) 



\^f) 



*J 



where '""^ 



are two orthonormal bases for the sub- 



fr) and \tpf 

systems A and B. Since the matrix occurs during Schmidt 
decomposition, we refer to it as the Schmidt matrix. The sin- 
gular value decomposition (S YD) of M provides the Schmidt 
coefficients. The eigenvalues {Ai} of the reduced density ma- 
trix are the squares of the singular values of M. It is very 
convenient, therefore, to express perturbative considerations 
on the ES, {^i = — log A^}, in terms of the M matrix. 
In the main text, we mentioned the following observation: 



When calculating IVI perturbatively, we find 
through examination of several examples that a 
contribution to M at some order can lead to a 
new ES level only if the contribution is not ap- 
pearing on the same row or the same column as a 
previous contribution which led to a new ES level 
at lower order. In physical terms, this means that 
we do not get new ES levels by applying a per- 
turbation only to the A block but need to apply 
perturbation terms to both blocks in a linked way. 



This observation is related to how the rank of the M ma- 
trix grows when the matrix M is being built up perturbatively. 
These somewhat abstract statements will be made concrete in 
these Supplements, by explicitly writing out the M matrix for 
several cases. These demonstrations of Schmidt matrices will 
also clarify the physical significance of this feature of per- 
turbative ES calculations, in terms of distinguishing between 
'bulk' and 'boundary-connected' contributions to the M ma- 
trix. 



Section will present the perturbative structure of the M ma- 
trix in a very general form, demonstrating the above statement 
(and additional features of M) under relatively mild assump- 
tions, without reference to an explicit Hamiltonian. This anal- 
ysis will also lead to a general 'proof of boundary-locality 
of the ES, for low perturbation orders. This complements 
the material in the main text, where our physical results on 
boundary-locality have been illustrated mainly through exam- 
ples on explicit physical systems and not generically. 

Section will present the perturbative calculation of the ES 
for the case of the Heisenberg ladder in some detail, for first 
and second order. This complements the main text where the 
main example was the XXZ chain, and where we focused on 
physical concepts at the expense of detailed formulas. In par- 
ticular, we will show for this explicit model how the M matrix 
grows perturbatively, leading to more and more ES levels at 
increasing orders. 



cialize to systems where KfRf^^ is hermitian.) Again, in- 
cluding these contributions would be straightforward but cum- 
bersome, and does not change any of the general conclu- 
sions below. For simplicity we will also assume that apply- 
ing KlKf_^^ on the unperturbed ground state any number of 
times generates product states which are also eigenstates of 
the unperturbed Hamiltonian ^^ "H^. 

The unperturbed product state |0) = Yli |0i) has the simple 
Schmidt decomposition |0) = |0)a(8'|0)b- Therefore, at order 
zero in perturbation theory, the M matrix has only the single 
element an = 1. 

First order — For first order perturbation theory, we can 
choose as basis for the A subsystem the states 



\0)a 
\0')a 



i,i + l e A 



PERTURBATIVE STRUCTURE OF ES AND OF M MATRIX: 
GENERAL FEATURES 

In this section we present in very general form the perturba- 
tive structure of the ES and of the M matrix, without relying 
on a particular example system and using only mild assump- 
tions about the form of the Hamiltonian. This also leads to a 
relatively general proof of boundary-locality, for low pertur- 
bation orders. We do not attempt or claim completeness of 
proof. Our purpose is to show why boundary-locality is ex- 
pected to hold as a general feature of the ES in gapped phases 
connected to a product state. 

We consider the system Hamiltonian to be of the form 



^ = E^» = E^° + ^E^^ 



(4) 



where e is a small parameter. The interaction is a sum of local 
terms, which we assume to be finite range but not necessarily 
single-site. The hamiltonian is translationally invariant (Hi = 
T-Li+i Vi), except at the system edges. 

The unperturbed ground state is taken to be a product state, 
as is true in all our examples. The unperturbed Hamiltonian 
(Si ^i*) ™^y ^^ either single-site or finite-range — in our 
example systems, T-Q are single-site (single-rung) terms in 
the Bose-Hubbard chain and the Heisenberg ladder, but of 
nearest-neighbor form in the XXZ chain. 

We specialize to the generic case of a nearest neighbor in- 
teraction of the form %[ — KjKf^-^ + h.c, where Kj, Kf 
are arbitrary operators acting on site i. This includes a large 
class of models, including the Bose-Hubbard chain and the 
XXZ chain. It is straightforward but cumbersome to gener- 
alize our results below to the case of next-neighrest-neighbor 
(or generic finite range) interactions, or to other local Hamil- 
tonians. 

For simplicity, we write below only contributions arising 
from the KlKf_^_^ term, and omit the terms or elements aris- 
ing from the hermitian conjugate. (Or, equivalently, we spe- 



and a corresponding basis set for the B block. There are 
as many 'bulk' states KlKf^-j^\0)A in this basis as there are 
bonds in A. (ci is the normalization.) The two "boundary 

states" 10') 



MB, 



are defined such that 



KJK]+^\0), 



with 



>s-c2|0')a®|0% 
j e A, J + le B 



i.e., \0')a{B) ^e the A(B) states produced when the perturba- 
tion H'j = K^K^j_^ acts across the boundary. We thus have a 
total of 1 + (i/i(B) ^ 1) + 1 basis states, where i/i(B) — 1 is 
the number of bonds in block A(B). 

[Our omission of the hermitian conjugate has simplified the 
number of "boundary terms" at first order. For a generic non- 
factorizable interaction of the form Kj j^i we would have 

^w+i|Oi)®|o,+i) = E™,„ 

is a basis for the local Hilbert space at site j.] 

With the above-described (1 + [L^^^-j — 1] + l)-dimensional 
bases for A(B) blocks, the M matrix at first order is of {Lb + 
1) X {La + 1) form 



la") ® la^^+i) where aj" 



(an 


ae 


ae ■ 


• ae 


\ 


ae 
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■ 





ae 





• 


• 






ae 





• 


• 





\o 





• 


• 


asse/ 



This matrix can be simplified by choosing the more conve- 
nient basis for the A block 



|0). 



ci 



Vl, 



1 ^ 

i.i+leA 



KlKl,\ 



10') 



and the corresponding set for block B. The idea here is that 
the Kf Kf^^\0) A states can be combined together as they all 
provide the same contribution. Going to this new basis first for 
the A block and then for the B block, we get the 'shrinking' 



/an 


ae 


ae ■ ■ 


• 


ae 





•• 


• 


ae 





•• 


• 



V 



033 £/ 



/an ai2e \ 
ae 
ae 



Vo 



033 £/ 



an ai2e 
aaie 
0336; 



so that the M matrix can be written as 




M 



with ai2 = ayJ{LA — 1) and a2i = ayJ{LB — 1). This 
"shrinking" property holds in general: if M has an isolated 
block with only one nonzero row or column, since its rank is 
trivially one it is possible to shrink the block to form a column 
or a row with only one nonzero element. 

The matrix M written above has a 2 x 2 block with the two 
Schmidt (singular) values: 



an 



2an 



-^21 2 



012021 2 
On 



(5) 
(6) 



In a first-order calculation, O(e^) contributions to the 
Schmidt values are not meaningful. The second singular value 
above is thus 'fake', and the 2x2 block has rank 1, up to 
meaningful order. The 0{e'^) correction to the first singular 
value is also meaningless. 

Thus, the ai2,02i appearing on the first row and column 
respectively do not give rise to any first order Schmidt value 
since they contribute only at order O(e^). This illustrates in a 
general setting, not referring to a particular model, that matrix 
elements appearing on the same row or column as previous 
contributions which gave lower order ES levels do not give 
new levels. Such terms only contribute to a renormalization 
of \Q) A{B)- We note that these non-contributing terms origi- 
nate from perturbative processes acting only on one of the two 
subsystems. 

In contrast, the element a33e forms a separate block by it- 
self, and thus contributes a new 0{e) ES level. This is an 
example of an element appearing perturbatively in a new row 
and new column of the M matrix and giving rise to a new 
ES level. This new ES level comes from the boundary terms 

|0')a(s). 

We thus have proved in a very general setting that, at first 
order, bulk terms disjoint from the boundary do not contribute 
to the ES and the only contribution comes from the boundary. 

Second order — At second order, we need to include ex- 
tra basis states. All or most of these new states are formed 
by applying T-L[ twice within the same block. This leads to a 
set of 0{e'^) contributions which in the M matrix fall on the 



r 



same row or column as an. As in first order, we can shrink 
these to a single new row and a single new column, and as in 
first order we can show that they do not provide any true ES 
contributions at this order Again, bulk contributions are seen 
not to contribute to the ES at low order. We omit such terms 
in the expressions written below. 

More interesting second-order contributions are those 
where T-L[ acts once in A and once in B. This can be writ- 
ten within the basis set used previously for first order The 
2x2 block of the M matrix, discussed previously at first or- 
der, takes the form 

'on 0126 
a2ie a22e^ 

The singular values from this 2x2 block are 



On 



''12 



2an 



-^21 2 



012021 



022O11 2 



On 



(7) 
(8) 



The extra value (as well as the renormalization to an = 1) 
is now relevant since it is of same order as the perturbation 
order. We thus can have a new contribution, but only if a22 7^ 
012021. 

This condition again reflects boundary locality. For second- 
order terms where the two applications i/' occur away from 
the boundary in A and B, and therefore away from each other, 
the contributions 'factorize' into two first-order applications 
of _ff' (encoded in the ai2, a2i terms). The physical processes 
that might cause a22 7^ 012021 are therefore appUcations of 
the perturbation connected to the boundary. (This intuitive 
argument is difficult to prove in full generality, but we have 
found it to hold in all the systems we have treated.) In our 
example systems, we also find that such a second-order con- 
tribution also requires the unperturbed Hamiltonian (T-Q) to be 
of longer range than single-site or single-rung terms. Thus the 
XXZ chain has such second order ES levels, while the Heisen- 
berg ladder and the Bose-Hubbard chain lacks such second 
order levels, as can be seen from Figures 1 and 3 in the main 
text. 

There are additional possibilities at second order, but they 
result in minor variations of the points discussed above. 

Thus, also at second order in perturbation theory, contri- 
butions to the ES are seen very generally to come only from 
boundary-connected processes. 



HEISENBERG SPIN LADDER 

In this section, we describe perturbative calculations of the 
ES in the spin-^ ladder with SU{2) Heisenberg interactions. 
We have performed calculations up to third order, but the 
expressions and notation become extremely cumbersome to 
write out, so we give details up to second order All the phys- 
ical points we want to demonstrate are already visible by this 
stage of the perturbative calculation. 



(1) 

(2) 



B 



L,-l ^A ;^A+1 ^A+2 



'l-g J„. 



FIG. 4. Spin ladder. The boundary between A and B regions is 
shown with a dashed vertical line. The rungs in block A and B are 
numbered from 1 to La and from La + 1 to La + Lb, respectively. 



Hamiltonian 



The Hamiltonian for the two-leg Heisenberg ladder is 



H = Hr 



ung 



H, 



•Jrung / ^ >Ji J, 



leg 

(1) o(2) 



•^'eg 2^ [ S^ ■ Dj^i 



<;;.(2) o(2) 



(9) 



where Jmng is the coupling between spins on the same rung, 
and J/eg is the coupling along each leg (Figure 4). The sub- 
script i indexes rungs, and (1),(2) refer to the two legs. Open 
boundary conditions are used, so that there is a single bound- 
ary between the blocks, shown with a dashed line in Figure 4. 
The number of rungs is La in block A and Lb in block B; 
total L = La + Lb- 

We are interested in the limit Ji^g <C Jmng where Hi^g can 
be treated as a perturbation. We will use Jieg = 1. The small 
parameter for perturbation is 



^leg I 'Jr 



XjJr- 



The eigenstates of the rung Hamiltonian Hmng are conve- 
niently expressed in terms of the single rung states. These 
are the singlet \s) and the three triplets |<"): 



1=^) 



1 



V2 



(I n> - 



in) 



it"> 



1 



V2 



(I n> + 1 ;t» 



The rung energy is 



1 / 

A ^1 



At 

4 '-'rung 



for the singlet and 



^ for the triplets. The difference is 1/e. 
The zeroth order ground state is the product state 



4 'J rung 



^ 



(0) _ 



ss ■ ■ ■ SS) — \S)i 



\S)2-- 



\S)L 



i.e., each rung is a singlet state. Other unperturbed eigenstates 
are those in which some of the rungs have been converted to 

r). 

To consider the result of applying Hieg on any state, it is 



convenient to know how Hi^g acts on a pair of rungs: 



m,g\t°t^) = -{\ss) + \t+t-) + \t-t+)) 
Hl,g\t+t-) = li-\t+t-) + \t't°)-\ss)) 

m,g\t-t+) = ^i-\t-t+) + \t't")-\ss)) 

Hieg\t"s) = ^\sn a = +,-,0 

i?ies|i^O = ^|tV) /? = +,- 

In the first row the sum is over a = +, 0, — and the sign 
function (" is (" ^ —1 if a = + , — and C" = 1 for a — 0. 



First order 

A single application of the perturbation operator creates a 
|i" ) |P ) pair in the sea of |s ) rungs. The wavefunction at first 
order is 



V,W=^(o)+ Y,c.o.t 



(0) 
ia 



with 



,,(0) 



2(l/e) 



-Ce/4 



Here tplJ is the state where the singlet pair i, i + 1 has been 
promoted to a |i"i" ) pair The unperturbed ground state -0'^*'^ 
appears without a subscript. 



V' 



(1) _ 



E(i^)i---k"*")''---i->L) 

i 

E(i»>i---l'*'">.---i»>i 

'-in 



s I--- t-r 



\S)L 



We introduce the following concise notation for states: all 
sites which are in \s) states are omitted and only the |i")'s 
are explicitly written, with the subscript indicating the site or 
the bond. Thus 



V'(^) = V' 



(0) 



^^r|i"t^). 



(10) 



Note that the norm deviates from unity by some O(e^) value, 
i.e., is correct at first order. 

We will now construct the Schmidt (M) matrix at first or- 
der The block spin component S% is a good quantum number; 
therefore we treat each spin sector separately. 



S^ = 1 sector ■ 



We start with the S*^ 



1 sector. The 



relevant basis for the two subsystems A and B is 



BA = {\t+)LA: 



Be 



{\t-)L^ + l}- 



(11) 



The triplets are on rungs neighboring the boundary between 
the blocks. The \s) rungs have been omitted. 

The Schmidt matrix in this sector has a single element, 
e/4, which gives the reduced density matrix eigenvalue A' = 
e^/16, and the entanglement level 



e 



■\og[{e/ Af 



-21og(e/4) 



The same result is obtained in the S*^ 



-1 sector. 



S\ = sector - 



The S\ — Q sector is more interesting. 



We choose the block bases 

where \s) denotes the unperturbed state of that subsystem, and 
we defined 



1 



»A/i 



M, 



A/B 



^c"|i"i"; 



(12) 



where Na/b = ^/^iLA/B — 1) is the norm. As in Section 
we are combining states from boundary-disconnected excita- 
tions into a single state, to 'shrink' rows and columns of the 
M matrix. With this choice the M matrix in this sector has the 
form 



1 


-fA/'A 





fA/B 














— ■! 



The isolated element —e/4 gives the same contribution to the 
ES as the S'^ = ±1 sectors, completing the S = 1 multiple! 
at ^ = —2 log(e/4). (Lowest triplet in Figure 3a of the main 
text). 

The remaining 2x2 block has two singular values, but 
one of them is of higher than linear order in e, and is thus 
a 'fake' singular value that does not contribute to this order of 
the perturbative calculation. In other words, the effective rank 
of the 2x2 block, at the order being considered, is 1 . 



Specifically, the reduced density matrix eigenvalues from 
this block are 

\o = l + ^AJ^l+Ul)e^ + 0{e^) (13) 

lb 



Ai 



256 



{NAMBYe^ + 0{e'') 



(14) 



The second (Ai) is 'fake', as explained above. 

Thus, we see here one example of the observation quoted 
in the Introduction: since the first-order contributions €Ma.b 
appear on the same row or same column as the zeroth order 
element 1, they do not lead to new ES levels. On the other 
hand, the isolated matrix element —e/4 appears on a new row 
and new column, and therefore is allowed to give rise to a new 
ES level, which it does. 

Note that the 'fake' eigenvalues are size-dependent: this 
reflects the physics that they arise from perturbations acting on 
the bulk and not those acting on an edge-connected location. 
To neglect terms like e^Ma-N'b in the thermodynamic limit, 
one has to use the correct order of limits. 

The size-dependent correction in Aq to its unperturbed 
value (unity) is also 'fake'. We will later see that this fake 
correction gets canceled at next order in the perturbation. 



2nd order 



The state up to second order can be written as 



V' 



(2) _ 



(- 


3 


^^{La 


+ Lb- 


-i))v^ 


,(0) 


e 
4 + 


^) 




' ee). 


e2 





16 



/•a I iCt lOi \ 



j^J^CC^l't^f's^t^thi (15) 

a.fS.i n 



We use the notation for total wavefunctions introduced in the 
previous subsection. In addition, here s" denotes a sequence 
of n rungs in state s. 

S\ = 1 sector — We start again with the 5^ = 1 sector, 
which at first order had the lone element e/4. We choose the 
following basis 

BA = {\t+)L.,\cj,'A}, SB = {|t-)L. + l,|0ij)}, 

where we defined 



Wa 



1 



a. 71 



(16) 



with Ma = V 1 + ji^A - 2), and analogously for the B 
block. When we omit a site subscript on a block ket, the sym- 
bols in the ket refer to the configurations of the rungs ending 
at the boundary. The rungs further from the boundary are of 
course understood to be singlets. 



10 



In this basis the Schmidt matrix is 



+ *) +^ 




from which the reduced density matrix eigenvalues are found 
to be 



^'^S+i^+^('+-^^+-^^^''+^^''^ ^^^^ 



A" = ^(ACiA/-s)^(e^) + 0(0. 



(18) 



We see that A' has acquired an 0{e^) correction, while the 
size-dependent 0{e'^) term is spurious. The A" is a com- 
pletely spurious eigenvalue; the meaningful rank of this 2x2 
block is 1, at this order. 

S*^ = sector — The S^ — sector itself has the block 
structure 



M, 
Mt 



where the triplet part Mt is the same as the matrix obtained 
above for the 5"^ = ±1 case. 

For the singlet part Mg, we use the A block basis 

BA = {\s),\cbA),\^l)} (19) 

where 



\<P1] 



^ \ a,i 



^^CC^\t'^t''s''th^ 



a,/3,i n 



withAf'J^ = ^3{La-2) + 1{La-2){La-3). We find 
the singlet part of the S^ = Schmidt matrix to be 



M.,= 



(l 



\ 



f5(AAl+AA| + 3) -{1 + \Wa 



-(! + \We 



-hS" 



\NaMi 



AJ^B 







\^l\ 



/ 



The term e^(A(4 + A/"^ + 3)/32 appears due to normalization. 
The reduced density matrix eigenvalues obtained from Ms 
are, for simplicity, written below only for the case La = Lb'- 



Ao = l 



16 



+ |(i^-l)e3 + 0(e^) 



(20) 



A+ = 



8192 ~ 



{La-1)\ -14 + 9La + 3L1 



±Aj2{LA~l){-iL\ + LA~Q)]e'' (21) 



The extensive 0{e^) term that appeared in Eq. (13) is no 
longer present in Ao, and the contribution 3e^/16 is genuine 
second order. The spurious eigenvalue Ai of Eq. (14) has been 
pushed from O(e^) to 0{e^) in the A± here. 



Perturbative expressions for ES levels 

We have elaborated on the perturbative calculations at the 
first and second orders. We do not show details of the next or- 
ders, but the perturbative expressions at the lowest few orders 
for the ES levels, and the leading orders of the corresponding 
entanglement eigenstates, are listed below: 

& = \ss) 



-21og(f) 
-21og(f^) 



-21og( 



1024" I 



rs)] 

\tH^)] 



\t-t+)), 



\t-t+))- 



where /3 = ± while a ~ ±,0. The notation for the eigen- 
vectors above is that only the states of the two rungs of the A 
block nearest the boundary are shown. 

A noteworthy feature is that there is no O(e^) contribution. 
This reflects an aspect of boundary-locality on which we have 
commented in detail in Section . 

The perturbative expressions for ES levels are shown as red 
straight lines in Figure 5, and are seen to describe the exact 
ES levels rather well at low orders. 




FIG. 5. Ladder entanglement spectrum, obtained using DMRG for 
a 20-rung ladder {La = Lb ~ 10), as a function of Jmng = £~^. 
The perturbative expressions are shown as red straight lines. 



